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^ ■ This talk presents some progress achieved in collaboration with A.Linde and 

D.LindeB! towards understanding the true nature of the global spatial structure 
of the Universe as well as the most general stationary characteristics of its time- 
dependent state with eternally growing total volume. 



In our opinion, the simplest and, simultaneously, the most general version of inflationary 
cosmology is the chaotic inflation scenario!. It can be realized in all models were the other 
versions™ of inflationary theory can be realized. Several years ago it was realized that 
inflation in these theories has a very interesting propertyBi which will be discussed in this 
talk. If the Universe contains at least one inflationary domain of a size of horizon (/i-region) 
with a sufficiently large and homogeneous scalar field 0, then this domain will permanently 
produce new /i-regions of a similar type. During this process the total physical volume of 
the inflationary Universe (which is proportional to the total number of /i-regions) will grow 
indefinitely. 

Fortunately, some kind of stationarity may exist in many models of inflationary Universe 
due to the process of the Universe self-reproduction!. The properties of inflationary domains 

1 Talk, presented at the TEXAS/PASCOS conference, Berkeley, December 1992. 
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formed during the process of the self-reproduction of the Universe do not depend on the 
moment of time at which each such domain is formed; they depend only on the value of the 
scalar fields inside each domain, on the average density of matter in this domain and on the 
physical length scale. This kind of stationarity, as opposed to the "stationary distribution" 
P c ((f),t) (the probability density to find the field inside a given Hubble domain at time t) 
cannot be described in the minisuperspace approach. 

In order to describe the structure of the inflationary Universe beyond one /i-region (min- 
isuperspace) approach one has to investigate the probability distribution P p ((f>,t), which 
takes into account the inhomogeneous exponential growth of the volume of domains filled 
by field Solutions for this probability distribution were first obtained inHS for the case of 
chaotic inflation. It was shown that if the initial value of the scalar field is greater than 
some critical value 0*, then the probability distribution P p (<j), t) permanently moves to larger 
and larger fields 0, until it reaches the field P , at which the effective potential of the field 
becomes of the order of Planck density (we will assume M p = 1 hereafter), where the 
standard methods of Quantum field theory in a curved classical space are no longer valid. By 
the methods used iniHit was impossible to check whether P p (<f>, t) asymptotically approaches 
any stationary regime in the classical domain < <\> v . 

Several important steps towards the solution of this problem were made by Aryal and 
Vilenkini, Nambu and SasakiS and Mijic0. Their papers contain many beautiful results and 
insights, and we will use many results obtained by these authors. However, the generalization 
of the results oM to non-trivial case of chaotic inflation appeared to be not an easy problemBH. 
Mijic0 did not have a purpose to obtain a complete expression for the stationary distribution 
PJ&,i). The corresponding expressions were obtained for various types of potentials V{(f) 
in0. Unfortunately, according to0, the stationary distribution P p ((Lt) is almost entirely 
concentrated at 3> P , i.e. at V(<f>) ^> 1, where the methods used inll3 are inapplicable. 

In this talk we will argue that the process of the self-reproduction of inflationary Universe 
effectively kills itself at densities approaching the Planck density. This leads to the existence 
of a stationary probability distribution P p {4>, t) concentrated entirely at sub-Planckian den- 
sities V{4>) < 1 in a wide class of theories leading to chaotic inflation. 

Let us consider the simplest model of chaotic inflation based on the theory of a scalar 
field minimally coupled to gravity, with the Lagrangian 

1 -R + ld^oy- v{4>) . (i) 
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Here G = M~ 2 = 1 is the gravitational constant, R is the curvature scalar, V(<p) is the 
effective potential of the scalar field. If the classical field is sufficiently homogeneous in 
some domain of the Universe (see below), then its behavior inside this domain is governed 
by the equations 

+ 3H<j) = -dV/d<p } (2) 
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Here H = a/a, a(t) is the scale factor of the Universe, k = +1, —1, or for a closed, open or 
flat Universe, respectively. 

The most important fact for inflationary scenario is that for most potentials V(0) (e.g., 
in all power-law V(<p) = g n <p n /n and exponential V(<p) = ge a( ^ potentials) there is an in- 
termediate asymptotic regime of slow rolling of the field and quasi-exponential expansion 
of the Universe. This expansion (inflation) ends at ~ e where the slow-rolling regime 
<C 3if (0)0 breakes down. 

If, as it is usually assumed, the classical description of the Universe becomes possible only 
when the energy-momentum tensor of matter becomes smaller than 1, then at this moment 
<9 M 0<9 M < 1 and V((p) < 1. Therefore, the only constraint on the initial amplitude of the 
field is given by V(<p) < 1. This gives a typical initial value of the field in the theory ([!]): 



where P corresponds to the Planck energy density, V(4> p ) = 1. 

During the inflation all the inhomogeneities are stretched away and, if the evolution of 
the Universe were governed solely by classical equations of motion (H), @, we would end up 
with extremely smooth geometry of the spatial section of the Universe with no primordial 
fluctuations to initiate the growth of galaxies and large-scale structure. Fortunately, the 
same Jeans instability which causes the growth of galaxies during the Hot Big Bang era 
leads to the existence of the growing modes of vacuum fluctuations during the inflation. 
The wavelengths of all vacuum fluctuations of the scalar field grow exponentially in the 
expanding Universe. When the wavelength of any particular fluctuation becomes greater 
than H~ x , this fluctuation stops oscillating, and its amplitude freezes at some nonzero value 
5cj)(x) because of the large friction term 3if0 in the equation of motion of the field 0. The 
amplitude of this fluctuation then remains almost unchanged for a very long time, whereas 
its wavelength grows exponentially. Therefore, the appearance of such a frozen fluctuation 
is equivalent to the appearance of a classical field <50(x) that does not vanish after averaging 
over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the cre- 
ation of more and more perturbations of the classical field with wavelengths greater than 
H~ x . The average amplitude of such perturbations generated during a time interval H~~ l (in 
which the Universe expands by a factor of e) is given by 



If the field is massless (or, better to say, as long as m 2 (0) <C if 2 (0)), the amplitude of 
each frozen wave does not change in time at all. On the other hand, phases of each waves 
are random. Therefore, the sum of all waves at a given point fluctuates and experiences 
Brownian jumps in all directions in the space of fields. 

The standard way of description of the stochastic behavior of the inflaton field during 
the slow-rolling stage is to coarse-grain it over /i-regions and consider the effective equation 



O ~ 0, 



P ' 



(4) 




(5) 
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of motion of the long-wavelength 

Here £(t) is the effective white noise generated by quantum fluctuations, which leads to the 
Brownian motion of the classical field (f>. 

This Langevin equation corresponds to the following Fokker-Planck equation for the 
probability density P c (<j),t) to find the field in a given point (which now means /i-region) 
after time t: 

at dt{ si 2 d<i> y [9, C ) + 3H(4>) C )' ( > 

The formal stationary solution (dP c /dt = 0) of equation (|7|) would be§ 

/ 3 \ 

P c ~ exp 



W(<f>)) ' 

Note, that the stationary solution is equal to the square of the Hartle-Hawking wave 
function of the Universelll. At first glance, this result is a direct confirmation of the Hartle- 
Hawking prescription for the wave function of the Universe. 

However, in all realistic cosmological theories, in which V(<f>) = at its minimum, the 
distribution (H) is not normalizable. The source of this difficulty can be easily understood: 
any stationary distribution may exist only due to a compensation of a classical flow of the 
field cf) downwards to the minimum of V(4>) by the diffusion motion upwards. However, 
diffusion of the field <fi discussed above exists only during inflation, i.e. only for <fi > 1, 
> ~ m2 ~ 10 -12 for m ~ 10~ 6 . Therefore @ would correctly describe the 

stationary distribution P c (4>) in the inflationary Universe only if V{4>) > 10 -12 ~ 10 80 GeV 
in the absolute minimum of V(0), which is, of course, absolutely unrealistic!. 

It was shown! that the solutions of this equation with the effect of "end of inflation" 
boundary properly taken into account are non-stationary (decaying). It is of no surprize 
because we didn't take into account yet the complicated and inhomogeneous expansion of 
the whole Universe with multiplicating number of /i-regions. The close view on the process 
of the growth of the volume of inflationary Universe reveals the strong resemblance with 
branching diffusion processes, with the role of branching particles being played by /i-regions 
and the the diffusing parameter being associated with the inflaton field inside an ^-region. 
If one wishes to describe the distribution of the inflaton field in the whole Universe and not 
only in one Hubble domain, one has to modify the Fokker-Planck equation by introducing 
the term corresponding to creation of more and more new "particles" during the diffusion 
process 0,100. 
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The mathematical model describing such behavior is the recently developed theory of 
a new type of branching diffusion processes^. Eq. (||) is, from this perspective, the for- 
ward Kolmogorov equation for the first moment of the generating functional of number of 
branching particles. 

There are two main sets of questions which may be asked concerning such processes. 
First of all, one may be interested in the probability P p (4>, t) to find a given field at a given 
time t under the condition that initial value of the field was equal to some 0(t = 0) = </>q. In 
what follows we will denote 0o as X- On the other hand, one may wish to know, what is the 
probability P p {Xi t) that the given final value of the field (or the state with a given final 
density p) appeared as a process of diffusion and branching of a domain containing some 
field x- Or, more generally, what are the typical history of branching Brownian trajectories 
which end up at a hypersurface of a given <p (or a given p)? 

It happens that if we do not take the Planck boundary seriously, this equation also 
doesn't have a true stationary solution (the solutions found by Nambu et. al. are heavily 
concentrated at the super-Planckian densities which is just another way to state that there is 
no real stationarity). This conclusion, however, may fail if we will treat the Planck boundary 
more carefully. There are different reasons to do this: 

1. Eq. (H) may have a slightly different form, which corresponds to the difference between 
stochastic approaches of Ito and Stratonovich. This difference is not important at densities 
much smaller than 1, but at V(<p) > 1 it may become significant. (Note, however, that this 
is not an unsolvable problem. One may just investigate modified equations as well. Two 
other problems are more fundamental.) 

2. Diffusion equations were derived in the semiclassical approximation which breaks 
down near the Planck energy density p p ~ = 1. 

3. Interpretation of the processes described by these equations is based on the notion 
of classical fields in a classical space-time, which is not applicable at densities larger than 1 
because of large fluctuations of metric at such densities. In particular, our interpretation of 
P c and P p as of probabilities to find classical field in a given point at a given time does 
not make much sense at p > 1. 

There is also another, more general, reason to expect the existence of the stationary 
solutions: as we will argue now, inflation kills itself as the density approaches the Planck 
density p p ~ M^. 

In our previous investigation we assumed that the vacuum energy density is given by V{<p) 
and the energy-momentum tensor is given by V((p) g pu . However, quantum fluctuations of 
the scalar field give the contribution to the average value of the energy momentum tensor, 
which does not depend on mass (for m 2 <C H 2 ) and is given byEl 

<T, U >= y^g»„ = \v 2 9flv . (10) 

The origin of this contribution is obvious. Quantum fluctuations of the scalar field <fi freeze 
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out with the amplitude j- and the wavelength ~ H 1 . Thus, they lead to the gradient 
energy density (d p 5(f>) 2 ~ H 4 . 

An interesting property of eq. fllOD is that the average value of the energy-momentum 
tensor of quantum fluctuations does not look as an energy-momentum tensor corresponding 
to the gradients of a sinusoidal wave. It looks rather as a renormalization of the vacuum 
energy-momentum tensor (it is proportional to g^). This means, in particular, that after 
averaging over all possible outcomes of the process of generation of long-wave perturbations, 
the result (for < T^ u >) does not depend on the choice of the coordinate system. However, 
if we are interested in local events (i.e. we are averaging only over short wavelengths), we 
will see long-wave inhomogeneities produced by the "frozen" fluctuations of the scalar field. 
This does not lead to any interesting effects at V <C 1 ((f) <C (ftp) since in this case V 2 <C V. 

However, at the density larger than the Planck density the situation becomes much more 
complicated. At V > 1 the gradient energy density ~ V 2 becomes larger than the potential 
energy density V((p). A typical wavelength of perturbations giving the main contribution 
to the gradient energy is given by the size of the horizon, I ~ if -1 . This means that the 
inflationary Universe at the Planck density becomes divided into many domains of the size 
of the horizon, density contrast between each of these domains being of the order of one. 
These domains evolve as separate mini-Universes with the energy density dominated not 
by the potential energy density but by the energy density of gradients of the field 0. Such 
domains drop out from the process of exponential expansion. Some of them may re-enter 
this process later, but many of them will collapse within the typical time H -1 . 

This can be effectively described by imposing some kind of absorbing (or reflecting, or 
elastic screen type) boundary conditions on our process of branching diffusion at some 
which should be close to the Planck boundary. Without going into details of derivation, we 
just announce here the resultsi. 

Our calculations have shown! that if our boundary conditions do not permit the field 
penetrate deeply into the domain > (f> p , the final results of our investigation do not depend 
on the type of the boundary conditions imposed (whether they are absorbing, reflecting, 
etc.). They depend only on the value of the field <fib where the boundary condition is to be 
imposed, and this dependence is rather trivial. We have argued that <pb ~ <fi p - Therefore 
we will assume now that the function P p ((f), t\x) satisfies boundary conditions corresponding 
to disappearing particles at <fi p treating P now just as a phenomenological parameter not 
necessarily corresponding to V(4> p ) = 1. In this case eq. (§) and its conjugate (backward 
Kolmogorov equation) are endowed by the following condition: 

P p ((P p ,t\x)=0. (11) 

One may try to obtain solutions of equations @ and its conjugate in a form of the fol- 
lowing series of biorthonormal system of eigenfunctions of the pair of adjoint linear operators 



6 



(defined by the equations below): 



P p (0 P ,t|x) = £e A ^( X R(0) ■ (12) 
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Indeed, this gives us a solution of eq. ([|) if 
The orthonormality condition reads 

> P 

V'aVX) Ti(X) rf X = 5 sj ( 15 ) 

In our case (with regular boundary conditions) one can easily show that the spectrum of 
Xj is discrete and bounded from above. Therefore the asymptotic solution for P p (<f>,t\x) (in 
the limit t — ► oo) is given by 

PpOMIx) = e Al Vi(x)Ti(0) • (l + 0(e- (Al - A2)t )) • (16) 

Here ^i(x) is the only positive eigenfunction of eq. fll3l) , Aj is the corresponding (real) 
eigenvalue, and 7i"i(0) (invariant density of branching diffusion) is the eigenfunction of the 
conjugate operator (|14j) with the same eigenvalue Ai. Note, that Ai is the largest eigenvalue, 
Re (Ai — A2) > 0. We found@ that in realistic theories of inflation the typical time of relax- 
ation to the asymptotic regime, t re i ~ (Ai — A2) - " 1 , is extremely small. It is only about a few 
thousands Planck times, i.e. about 10~ 40 sec. This means, that the normalized distribution 

P p (0,t| X ) = e- Alt P P (0 P ,t|x) (17) 
rapidly converges to the time-independent normalized distribution 

PMx) = P P (<P,t^ oo\x) = MX) 7Tl(0) • (18) 

It is this stationary distribution that we were looking for. The remaining problem is to 
find the functions V'i(x) an d (</>)> an d to check that all assumptions about the boundary 
conditions which we made on the way to eq. flIE|) are actually satisfied. 

After some calculations we came to the following expression for P p (<p\x) (note that the 
functions $(•) are essentially the same, the only difference is the argument): 

= ^ {wvm - ) * Wx)) Hzm (19) 
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Here the normalization constant C should be determined from ([151) . Using the WKB ap- 
proximation, we have calculated! $(z(</>)) for a wide class of potentials usually considered in 
the context of chaotic inflation, including potentials V ~ <p n and V ~ e Qfl \ 

The solution we have found features interesting properties. First of all, it is concentrated 
heavily at the highest allowed values of the inflaton field. Despite the promising exponential 
prefactor in (|T9|) which looks like a combination of the Hartle-HawkingB and tunneling)!! 
wavef unctions, the dependence of &(z((fi)) on appeared to be even more steep that those 
exponents. This function falls exponentially (with the rate governed by the Planckian ener- 
gies) towards the lower values of the inflaton field and strongly overwhelms the dependence 
of the exponent in front of it. Only near the "end of inflation" boundary the solution fllTf ) 
reveals something familiar — the dependence on the initial field x becomes similar to the 
square of tunneling wavefunction (simultaneously, the corresponding dependence on the final 
field cancels out). And, as we have already mentioned, the relaxation time is very small 
(which is also a consequence of the fact that the dynamics of the self-reproducing Universe 
is governed by the maximal possible energies). The stationary distribution found ini is not 
very sensitive to our assumptions concerning the concrete mechanism of suppression of pro- 
duction of inflationary domains with > <p p . We hope that these results may show us a way 
towards the complete quantum mechanical description of the stationary ground state of the 
Universe. 

A.M. is grateful to A.Linde and D.Linde for fruitful and enjoyable collaboration on this 
project and to A.Starobinsky, S.Molchanov and A.Vilenkin for discussions. 
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